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ABSTRACT 

The Rich-bum/Quick-mix/ Lean-bum (RQL) combustor has been 
identified as a potential gas turbine combustor concept to reduce NO x 
emissions in High Speed Civil Transport (HSCT) aircraft. To 
demonstrate reduced NO x levels, cylindrical flametube versions of RQL 
combustors are being tested at NASA Lewis Research Center. A critical 
technology needed for the RQL combustor is a method of quickly mixing 
by-pass combustion air with rich-bum gases. 

In this study, jet mixing in a cylindrical quick-mix section was 
numerically analyzed. The quick-mix configuration was five inches in 
diameter and employed twelve radial-inflow slots. The numerical 
analyses were performed with an advanced, validated 3-D Computa- 
tional Fluid Dynamics (CFD) code named REFLEQS. Parametric 
variation of jet-to-mainstream momentum flux ratio (J) and slot aspect 
ratio was investigated. Both non-reacting and reacting analyses were 
performed. 

Results showed mixing and NO x emissions to be highly sensitive 
to J and slot aspect ratio. Lowest NO x emissions occurred when the 
dilution jet penetrated to approximately mid-radius. The viability of 
using 3-D CFD analyses for optimizing jet mixing was demonstrated. 

NOMENCLATURE 

A Pre-exponential Factor 

D Diameter of Quick-Mix Section 

D/, Hydraulic Diameter 

El Emission Index 

E/R Activation Energy/Gas Constant 

J Jet-to-Mainstream Momentum Flux Ratio 

m, Mass Flow in Each Cell i 

m j Jet Mass Flow 

n\„ Mainstream Mass Flow 

or Mass Weighted Standard Deviation of Temperature 

n Optimum Number of Slots 

NO x Oxides of Nitrogen 

T fl pg Mass-Weighted Average Temperature 

Ti Temperature in Each Cell i 

u' RMS of U Velocity 

U Averaged Axial Velocity 

v' RMS of V Velocity 

V Averaged Radial Velocity 


INTRODUCTION 

In order to meet the growing need for faster transportation, High- 
Speed Civil Transport (HSCT) aircraft and associated propulsion 
systems have been under study in recent years. One major concern that 
has surfaced concerning HSCT engines is their impact on deteriorating 
the earth's ozone layer. Using current technology, a fleet of HSCT 
aircraft would produce large amounts of oxides of nitrogen (NO x ) while 
cruising in the stratosphere. Such high levels of NO X/ through a series 
of well known reactions, would drastically reduce ozone levels. In order 
to reduce NO x emissions, technology must be developed to design 
advanced, low emission combustors. 

One combustor concept that has been identified as a leading 
candidate to reduce NO x emissions is the Rich-burn/Quick-mix/ Lean- 
bum (RQL) combustor. Originally conceived and developed for 
industrial combustors (Mosier and Pierce, 1980 and Pierce et al. , 1980), 
the RQL concept utilizes staged burning, as shown in Figure 1. 
Combustion is initiated in a fuel rich zone at equivalence ratios between 
1.2 and 1.8, thereby reducing NO x formation by depleting the available 
oxygen. Bypass combustion air is introduced in a quick-mix section and 
lean combustion occurs downstream at an overall equivalence ratio 
between 0.5 and 0.7. The quick-mix section usually has a smaller 
geometric cross-section area than the rich bum zone in order to prevent 
backflow and enhance mixing. 



Fig. 1. Industrial Rich-bum/Quick-mix/ 

Lean-burn (RQL) Combustor (Pierce et al . , 1980) 

Perhaps the single most important issue of the RQL concept is the 
design of the quick-mix section. For previous laboratory combustors, 
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Tacina (1990) has shown RQL NO x levels to be higher than lean, 
premixed combustor NO x levels. The higher NO x emissions for RQL 
was attributed to stoichiometric burning in the quick-mix section, thus 
emphasizing the need for optimized rapid mix concepts. Indeed, 
Nguyen et al (1989) have shown that if instantaneous mixing is 
assumed in the quick-mix section, low NO x emission index can be 
obtained at HSCT cruise flight conditions. Hence, one challenge of the 
RQL concept is to identify quick-mix sections with rapid mixing. 

This study sought to investigate the influence of jet-to-mainstream 
momentum flux ratio (J) and slot aspect ratio (SAR) on mixing 
effectiveness in a RQL flametube combustor to be tested at NASA Lewis 
Research Center (LeRC). Conventionally, dilution air in can combustors 
has been introduced through radial inflow holes. According to 
Holdeman's correlation (Holdeman et al, 1987), optimum mixing occurs 
when the following expression is satisfied: 

O' 

C 

where 

n = optimum number of holes 

C = experimentally derived constant -2.5 

J = jet-to-mainstream momentum flux ratio (p jVj 2 / p«V„ 2 ). 

Unfortunately, this correlation was developed for circular holed 
dilution jet mixing and jet mass flow-to-mainstream mass flow ratios 
(my/rru) of approximately 0.5. The RQL combustor requires my/m* of 2.0, 
thus necessitating slots instead of holes around the can's perimeter. 
The design of slots for optimum mixing needs further investigation, and 
was the focus of this study. 

CFD CODE 

The approach in this study was to perform 3-D numerical computa- 
tions on a cylindrical quick-mix section. The goal of the study was to 
provide improved understanding of slot injection and mixing. An 
advanced CFD code, REFLF.QS, was used to perform the computations. 
REFLEQS was developed by CFD Research Corporation (Przekwas et 
al. , 1990 and Smith et al, 1988) to analyze turbulent, reacting flows. 
The basic capabilities/methodologies in REFLEQS include: 


After reviewing the time scales for heat release at cruise-type 
conditions in RQL mixers, it was determined reaction rates were much 
faster than mixing rates. Hence, the combustion process was considered 
mixing controlled, and instantaneous reaction rates for heat release 
were assumed. 

When rich burn gases (composed of equilibrium concentration of 
CO, CO 2 , H 2 , H 2 O and N 2 ) were mixed with air, they were assumed to 
react according to the equation: 

0.0570 CO + 0.0197 CO 2 + 0.0383 H 2 + 0.0680 H 2 0 + 0.1371 N 2 
+ 0.1047 02 + 0.3953 N 2 

-> 0.0767 CO 2 + 0.1063 H 2 0 + 0.7123 N 2 + 0.0565 0 2 ( 2) 

Accordingly, any CO concentration that remained in the mixer 
exit was the result of unmixedness. 

NO x MODEL 

It was assumed that the NO x reactions did not contribute to the 
overall heat release in the combustor, thus allowing the NO x reactions 
to be "decoupled" from the heat release reactions. NO x was calculated 
as a passive scalar after the computation of the reacting flowficld. 

A simple Zeldovich reaction scheme was used to model the NO x 
formation. According to the mechanism, NO formation can be described 
by: 

N 2 + O ++ NO + N (3) 

and 

O 2 + N++NO + O. (4) 

The first reaction is much slower than the second one and hence 
controls the rate of NO formation. If the concentration of NO is much 
smaller than the corresponding equilibrium value, the rate equation for 
NO can be written as: 


1. solution of two and three-dimensional Navier- 
Stokes equations for incompressible and compres- 
sible flows; 

2. cartesian, polar, and non-orthogonal body-fitted 
coordinates; 

3. porosity-resistivity technique for flows with 
internal blockages; 

4. fully implicit and strongly conservative for- 
mulation; 

5. three differencing schemes: upwind, hybrid, 
and central differencing with damping terms; 

6. standard (Launder and Spalding, 1974) and 
extended (Chen and Kim, 1987) k-z turbulence 
models, the two-scale turbulence model of Kim 
and Chen (1988), and the low-Reynolds number 
k-z model of Chien (1985); 

7. instantaneous, one-step, and two-step combustion 
models; 

8. modified form of Stone's Strongly Implicit 
Solver; and 

9. pressure-based solution algorithms including 
SIMPLE and a variant of SIMPLEC. 

REFLEQS has undergone a considerable amount of systematic 
quantitative validation for both incompressible and compressible 
flows. Over 30 validation cases have been performed to date, and good- 
to-excellent agreement between benchmark data and predictions has 
been shown (Smith et al., 1989; Ratcliff and Smith, 1989; and Avva et 
al . , 1990). The good agreement gives confidence in the numerical results 
of this study. 


WOI-jcwuioj (5) 

dt 

Approximating the concentrations of N 2 and O by the local equi- 
librium values, the rate equation is given by 


d (NO ) 
dt 


- A e 



(6) 


From Quan et al, (1972), the rate constants were determined to be: 


A = 5.74 x 10 4 .laoles 

(7) 

secy m 3 


£. = 6.7x10"°K 

(8) 


R 

The NO x model was calibrated against the experimental results 
of Anderson (1975), who used a premixed, prevaporized laboratory 
combustor. The REFLEQS test case consisted of premixed propane and 
air reacting in a straight channel and instantaneous heat release. The 
rate constants had to be modified to give good agreement with 
Anderson’s data. The final constants used in this study were: 

/\ = 3.3xio 4 -L .[mki (9> 

sec y 

£= 1.03x 10 l< < 10) 

R 
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Figure 2 shows the computed results compared to Anderson's data 
of Emission Index (El) as a function of adiabatic flame temperature. 
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Fig. 2. Calibration of NO x Model in REFLEQS 


upstream and downstream of the slot. The grid in the radial direction 
was non-uniform with greater density near the combustor wall (power 
expansion of 1.2). The grid in the transverse direction was uniform in 
the slot, and slowly expanding away from the slot. The baseline slot 
was represented by a 6x6 mesh. As will be discussed in grid indepen- 
dence studies, this rather coarse grid is not grid-independent, but it docs 
capture all of the relevant flow features. For comparative studies, it 
was felt sufficient. 

Numerical Details 

The numerical details of the calculations included: 

1. Whole field solution of u-momentum, v- 
momentum, w-momentum, pressure correction, 
turbulent kinetic energy k, dissipation rate e, 
total enthalpy, and mixture fraction; 

2. Upwind Differencing for parametric studies of J 
and central differencing for parametric studies 
of SAR; 

3. Variable Fluid Properties (i.e. temperature 
dependency of specific heat, laminar viscosity, 
etc.); 

4. Adiabatic Walls; 

5. Standard k-t Model with wall functions; 

6. Turbulent Prandtl number of 0.9; 

7. Instantaneous heat-release model; and 

8. Six active chemical species. 


DETAILS OF NUMERICAL MODEL 


A geometry compatible with the NASA LcRC flametube 
combustor was selected for analysis. The geometry, numerical grid, 
numerical details, boundary conditions, grid independence, and 
convergence criteria are discussed below. 


Gfomeby 

The geometry of the numerical model consisted of three com- 
ponents: an inlet pipe, converging section and a quick-mix section (see 
Figure 3). The inlet pipe was 0.152 m (6.0 in.) in diameter and 0.076 m 
(3.0 in.) in length. The inlet pipe converged into the quick-mix section 
which was 0.127 m (5.0 in.) in diameter (D). The length of the quick- 
mix section was 0.333 m (13.0 in.). In reality, the length of actual quick- 
mix section hardware is approximately 6 inches, but the computational 
domain is extended for better understanding of NO x formation and to 
eliminate flowfield contamination by exit boundary conditions. 



Fig. 3. Schematic of Quick-Mix Geometry 


Twelve slots were located symmetrically around the perimeter of 
the quick-mix section. The axial location of the slot centerline was 
0.076 m (3.0 in.) from the inlet of the quick-mix section. The baseline 
slots were rectangular in shape with an aspect ratio of four and aligned 
in the streamwise direction. Three variations in slot aspect ratio 
(SAR) were tested: 1, 4 and 16. 


Due to geometric symmetry, only one slot was modeled with 
planes of symmetry set up halfway between adjacent slots. This 
allowed greater grid resolution and reduction of computer turnaround 
time. The rO domain was reduced to a pie section with a central angle of 
thirty degrees. 



GRID: 



Fig. 4. Baseline Grid 


Boundary Conditions 


Grid 

A baseline grid of 9,216 cells (32x16x18 in x, r, 0 directions) was 
selected and used for modeling the mixer. The grid is shown in Figure 4. 
Note that the origin of the coordinate system is located at the center of 
the slot. The axial grid spacing is dense near the slot, and gets coarser 


Mainstream Boundary . At the mainstream inlet boundary, 
propane and air are assumed to have completely reacted at an 
equivalence ratio ($) of 1.6. The species and temperature of the 
reaction products were taken from the JANNAF-standard rocket code 
named One Dimensional Equilibrium (ODE) (Nickerson el al. , 1989). 
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Velocity and pressure were obtained from experimental test plans for 
the RQL flametube combustor. A uniform velocity profile was assumed 
with a turbulent intensity typical of primary zones in gas turbine 
combustors and a turbulent length scale corresponding to a turbulent 
viscosity 1000 times greater than laminar viscosity. The mainstream 
inlet conditions (at 6.0 in. diameter) were: 


Axial velocity 

Temperature 

Density 

Composition (mass fraction) 

Turbulent intensity (u'/U) 
Turbulent length scale (lj/DoJ 


= 355 m/s 
= 2221 °K 
= 232 kg/m 3 

= 0.134 CO, 0.068 CO 2 , 0.006 H 2 , 
0.096 H 2 0, 0.696 N 2 
= 50% 

= 0.02 


Since equilibrium NO x levels are very low for ^ of 1.6 (*= 4 ppm), no 
NO x was assumed in the mainstream inlet. 


jet Inlet . The composition at the dilution jet inlet was assumed to 
be air. A uniform velocity profile was assumed and turbulent properties 
were selected using the same logic as discussed for mainstream 
turbulence. The jet inlet flow conditions were: 


Mass flux ratio (my/mJ 
Jet temperature 
Density 

Composition (mass fraction) 
Turbulent intensity (v'/V) 
Turbulent length scale (l t /D/,) 


1.94 
81 1°K 
635 kg/m 3 
0.232 O 2 , 0.768 N 2 


10% 

0.13 


The momentum flux ratio (J) was varied parametrically from 16 to 
64 by variation of jet velocity from 120 m/s to 240 m/s. The jet velocity 
variation corresponded to liner pressure drops (AP/P) of 3 to 12 percent. 
For each jet velocity, the slot flow area was modified to maintain 
constant jet flow. 

Exit Boun dary. The exit boundary condition was a zero gradient 
boundary condition. 


S. The transverse boundaries were assumed 
to be symmetry planes. These boundaries were also tested for possible 
outflow by setting them to be periodic boundaries. No discernable 
difference was seen between cases with symmetric and periodic 
transverse boundaries. 


The combustor wall was treated as a no-slip 
adiabatic wall (zero enthalpy gradient). Wall functions were used for 
the calculations of wall shear stress and near-wall turbulent quantities 
( k and e ). 


Centerline . The computational boundary at the centerline was 
assumed to be a symmetry plane. 


Two different sizes of grids were run to test grid independence: 
9,216 and 52,650 cells. The finer grid was obtained by increasing the 
grid density by -75% in all three directions. Comparison of the two 
grids is shown in Figure 5. 

Computational results from the two grids are presented in Figure 6 
for a momentum flux ratio (J) of 32.0. The isotherms in an rx plane 
through the jet centerline are shown and compared. Qualitatively 
they exhibit similar features, although the jet penetrated a little 
further in the case of the fine grid. Isotherms are also shown for two 
axial planes: x/D = 0.0 and 2.0. The isotherms at x/D = 2.0 show 
slightly higher temperatures (~22°K) for the fine grid. Also, the cold 
region in the fine grid solution is located closer to the centerline, 
indicating greater penetration. However, overall the coarse grid 

solution is very similar to the fine grid solution, 

Based on this grid-independence study, it appears the coarse grid 
captures the overall physics of the problem, and can be used to 
qualitatively compare quick-mix designs. 



GRID: x = 32 9216 CELLS 

r = 16 
9 = 18 


Section A-A 




Section B-B 



Fig. 5. Comparison of Coarse and Fine Grids 
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Fig. 6. Comparison of Isotherms for 
Coarse and Fine Grids: J=32 


Convergence 

The summations of all error residuals were reduced five orders of 
magnitude, and continuity was conserved in each axial plane. 
Typically, convergence required approximately 150 iterations as shown 
in the Figure 7. The relaxation on the velocity components (u and v 
only) was continuously varied during the run through a user specified 
input file. The repeated variation of relaxation allowed resolution of 
different scales of numerical error. This was found to speed up conver- 
gence by a factor of six compared to constant relaxation. Approximately 
3 CPU hours were required on an Alliant FX/8 mini-supercomputer 
(operating on one computational clement). Fine grid calculations took 
approximately 500 iterations and 40 CPU hours. For comparison, the 
ALLIANT computer speeds are ~20 times slower than a CRAY X-MP. 

RESULTS 

Parametric numerical tests were performed for jet-to-mainstream 
momentum flux ratios (J), for both non-reacting and reacting gases. 
Parametric variation of slot aspect ratio (SAR) was also studied. 
Discussion of the findings are reported below. 


Variation of J; Non-Reacting JEio.iv 

Five jet-to-mainstream momentum flux ratios were parametri- 
cally tested: 16, 32, 40, 48, and 64. All other flow conditions were held 
constant, including mass flow ratio (jet-to-mainstream) at 1.94. To 
maintain a constant mass flow, the slot size was changed for each J. 
The slot aspect ratio was held constant at four, and was always 
centered at the same location. The same number of grid cells were used 
in all cases. However, since the slot size was changing, the grid density 
had to be slightly altered for each case. This variation is thought to 
have a minimal effect on the results discussed below. 

Computed temperature contour maps are presented at x/D of 1 .0, as 
shown in Figure 8. The radial location of the lowest temperatures 
indicates the penetration location of the cold jet. As expected, 
increased jet penetration can be seen for larger values of J. Best mixing 
seems to occur when the jet penetrates to approximately mid-radius. J 
of 32 and J of 40 appear to be optimum mixers. For comparison, the 
optimum J is 45.6 using Eq. (1). 
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Fig 8. Temperature Contour Maps for Non-Reacting 
Conditions: x/D=1.0 




Fig 10. Temperature Contour Maps for Reacting 
Conditions: x/D=1.0 
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For a more quantitative comparison of mixing effectiveness, the 
mass-weighted standard deviation of temperature (aj) was calculated 
for each case. Gj was defined as: 

Xm,(7'rT„,) ? (") 

Ira. 

<r r = 

Tmg 

In Figure 9, Oj is presented versus j. It can be seen that J of 32 has 
the lowest Gj at x/D = 2.0. Underpenetration is worse than over- 
penetration in terms of cj. 



Fig. 9. Mixing Effectiveness: Non Reacting Flow 
Variation of J: Reacting Flow 

The same cases were analyzed as discussed above, except chemical 
reaction was turned on. Due to reaction, the overall mass-averaged exit 
temperature increased from 1301°K for non-reacting flows to 1790°K for 
reacting flows. Figure 10 shows temperature contour maps for the 
reacting cases one diameter downstream of the jet center. From this 
figure, it appears that J = 40 is the best mixer. This can be further 
elucidated by looking at the mixing effectiveness (aj) shown in Figure 
11. Figure 11 shows J = 40 to be the best mixer at x/D =1.0 and 2.0. 



Fig. 11. Mixing Effectiveness: Reacting Flow 

In addition to mixing effectiveness, another important criterion 
for evaluation of quick-mix sections is combustion efficiency. In 
particular, CO concentrations should be essentially eliminated from 
the combustor exit. The CO emission level in various axial planes 
downstream of the dilution jet is displayed in Figure 12. For all cases 
except J = 16, it can be seen that the CO species has been oxidized (to 
CO 2 ) by x/D = 0.25. For J = 16, unreacted CO remains in the flowfield 
even at x/D = 2.0. This is due to jet underpenetration, thus allowing 
rich burn gases (containing CO) on the centerline to pass through the 
quick-mix section without contact with dilution air. 



Fig. 12. Predicted CO Emissions 

The NO x results are presented in terms of Emission Index (El) in 
Figure 13. For the optimum case (Jopt =40), El is 2.9 at x/D = 2.0. 
Significant increase in El is predicted as J is increased or decreased from 
the optimum value. For J greater than J opl , jet overpenetration causes jet 
backflow on the centerline, resulting in higher NO x emissions. For] less 
than J opt , underpenetration of the jet results in reaction (and high 
temperatures) on the combustor centerline. 
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levels are less than those for SAR of 4, but chemical reaction and NO x 
formation is still occurring at x/D of 1 due to jet underpenetration. This 
is evidenced by the steep slope of the NO x curve at x/D of 1. A similar 
effect of delayed NO x formation on the centerline caused by jet 
underpenctration was shown in Figure 14 for J of 16 and SAR of 4. Hence, 
the best SAR is 4, with jet overpenetration for SAR of 16 (and 
corresponding higher NO x levels). 

These results suggest the importance of SAR on NO x emissions. As 
was shown earlier for J variation, the jet must penetrate to ap- 
proximately mid-radius for optimum mixing, and hence lowest NO x . 



Fig. 13. Predicted NO x Emission Index 

Figure 14 shows NO x concentrations convected out of each axial 
plane. Except for J = 16, all the cases show very little NO x formation 
downstream of x/D = 1.0. This indicates that high temperature zones 
are no longer existent. For the J = 16 case, NO x formation is increasing 
significantly downstream of x/D = 1.0, indicating high temperature and 
chemical reaction is still taking place. 



Axial Distance (x/D) 

Fig. 14. History of NO x Formation in Mixer 


Variation of Slol Aspect Ratio 

Three slot aspect ratios (SAR) were numerically analyzed: 1, 4 
and 16. The long dimension was aligned in the mainstream flow 
direction. The numerical grid was slightly modified for each slot, and 
central differencing was employed for increased accuracy. The jet-to- 
mainstream momentum flux ratio and mass flow ratio was maintained 
constant at 32 and 1.94 respectively. 

Isotherms in the rx plane through the jet centerline are shown in 
Figure 15. As SAR increased, jet penetration increased (as seen in Figure 
15). This is due to reduced flow blockage as SAR is increased. Figure 16 
shows the effect of SAR on NO x emissions. For SAR of 1, predicted NO x 
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Fig. 15. Isotherms in rx Plane Through Jet Centerline: 
Effect of Slot Aspect Ratio (SAR) 


12 


N0 V Emission Index 



Fig. 16. Effect of Slot Aspect Ratio (SAR) 
on NO x Emissions 


CONCLUSIONS 

The overall conclusions of this study were: 

1. The viability of using 3-D CFD to model and 
screen quick-mix concepts of low emission 
combustors was successfully demonstrated. 

2. A five-inch diameter quick-mix section 
compatible with the NASA LeRC Low Emission 
Combustor Program was numerically analyzed. 

The configuration consisted of twelve, radial- 
inflow slots uniformly distributed around the 
perimeter of the quick-mix section. Optimum 
mixing for non-reacting flow occurred for a jet-to- 
mainstream momentum flux ratio (J) between 32 
and 40. For reacting flow, the NO x emission 
index was shown to be highly sensitive to J, 
with the lowest value of 2.9 calculated for J of 
40 (at x/D = 2.0). 

3. The numerical results suggest that slot aspect 
ratio has a pronounced effect on jet penetration 
and mixing effectiveness. Conventional 
correlations for optimum mixing effectiveness 
for holes may not be applicable for slots. 
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